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In years to come, we shall see substantial progress in our understanding of nuclear struc- 
ture - a rich and many-faceted field. An important element in this task will be to extend 
the study of nuclei into new domains. The journey to "the limits" of isospin, angular mo- 
mentum, and mass and charge is a quest for new and unexpected phenomena which await 
us in uncharted territories. What is extremely important from a theoretical point of view 
is that the new data are also expected to bring qualitatively new information about the 
effective nucleon-nucleon interaction and hence about the fundamental properties of the nu- 
cleonic many-body system. The main objective of this presentation is to discuss some of the 
, challenges and opportunities for nuclear structure research with radioactive nuclear beams. 

O ' 

^ \ §1. Introduction 

The range of unstable nuclei accessible with radioactive nuclear beam (RNB) 
facilities opens up enormous opportunities for the study of nuclear structure, nuclear 
astrophysics, and fundamental interaction physics EP'H^cP. Intriguing possibilities 
occur both at the drip lines and in the long iso-chains of nuclei between the valley 
of stability and the extremes of nuclear existence. Exotica in the latter region are 
almost sure to appear since the mean field in weakly bound neutron-rich nuclei is 
modified relative to nuclei near stability, since new types of correlations are likely 
to occur, and since reduced binding should modify residual interactions among the 
outermost nucleons. Between the regions of known and near-drip-line nuclei lies an 
extensive zone where studies will reveal much about the microscopy of the nucleus. 

In this paper, we intend to discuss - rather briefly - several theoretical chal- 
lenges related to the mean- field description of exotic nuclei. We focus on several 
recent developments and trends while the reader is referred to the reviewtf for more 
discussion on earlier studies and aspects not covered here. 

§2. Physics of Neutron-Rich Nuclei 

One of the main avenues addressed by radioactive ion beams is the evolution of 
nuclear structure as a function of neutron-to-proton asymmetry. 

Figure 1 shows various domains of nuclear matter, important in the context of the 
RNB program. The range of neutron excess, (N — Z) / A, in finite nuclei is from about 
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Fig. 1. Diagram illustrating the range of isoscalar densities and neutron excess, (N — Z)/A, impor- 
tant in various areas of nuclear physics and astrophysics (based on Ref. If). 



-0.2 (proton drip line) to 0.5 (neutron drip line). Of course, by using radioactive 
beams (and radioactive targets), one will be able to explore new regions of Fig. 1. In 
particular, the new-generation RNB facilities, such as the Rare Isotope Accelerator, 
will provide a unique capability for accessing the very asymmetric nuclear matter 
and for compressing neutron-rich matter approaching density regimes important for 
supernova and neutron star physics. 

From a theoretical point of view, exotic neutron-rich nuclei far from stability 
are of particular interest. They offer a unique test of those components of effec- 
tive interactions that depend on the isospin degrees of freedom. In many respects, 
weakly bound nuclei are much more difficult to treat theoretically than well-bound 
systems!!*. For weakly bound nuclei, the Fermi energy lies very close to zero, and 
the particle continuum must be taken into account explicitly. 

2.1. Uncertain Extrapolations 

Experimentation with radioactive nuclear beams is expected to expand the range 
of known nuclei. That is, by going to nuclei with extreme N/Z ratios, one can 
magnify the isospin-dependent terms of the effective interaction (which are small in 
"normal" nuclei). The hope is that after probing these terms at the limits of extreme 
isospin, we can later go back to the valley of stability and improve the description 
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Fig. 2. Predicted two-neutron separation energies for the even-even Sn isotopes using several 
microscopic models based on effective density-dependent nucleon-nucleon interactions (from 
Refs.rJ>>0>). 



of normal nuclei. In addition to nuclear structure interest, the understanding of 
effective interactions in the neutron-rich and proton-rich environment is important 
for astrophysics and cosmology. 

Figure 2 illustrates difficulties with making theoretical extrapolations into neu- 
tron-rich territory. It shows the two-neutron separation energies for the even-even 
Sn isotopes calculated in several microscopic models based on different effective 
interactions. Clearly, the differences between forces are greater in the region of 
"terra incognita" than in the region where masses are known (an effect half-jokingly 
called "asymptotic freedom"). Therefore, the uncertainty due to the largely unknown 
isospin dependence of the effective force (in both particle-hole and particle-particle 
channels) gives an appreciable theoretical "error bar" for the position of the drip 
line. Unfortunately, the results presented in Fig. 2 do not tell us much about which 
of the forces discussed should be preferred since one is dealing with dramatic extrap- 
olations far beyond the region known experimentally. However, a detailed analysis 
of the force dependence of results may give us valuable information on the relative 
importance of various force parameters. Moreover, it is obvious that the experi- 
mental data on nuclei near stability do not sufficiently constrain results for exotic 
nuclei. Therefore, further measurements in neutron-rich systems are of vital im- 
portance for the development of theoretical descriptions. We note in passing that 
attempts to derive the effective interactions from first principles may lead to even 
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larger uncertainties for exotic nuclei (see recent analysis in Ref.B*). 

Figure 3 demonstrates that one does not need to go all the way to the neutron 
drip line to see deviations from well-established trends. Indeed, the differences be- 
tween various models describing masses of Sn isotopes show up just above N=82. 
Interestingly, experimental data around 132 Sn exhibit an unusual pattern as one 
crosses iV=82. Namely, there is a striking asymmetry in the position of 2f levels in 
iV=80 and 82 isotopes of Sn and Te, and the B(E2; 0+ -> 2|) rate in 136 Te stays 
unexpectedly low if , defying common wisdom that the decrease in E + in open-shell 

nuclei must imply the increase in B(E2;0 + — ► 2+). So far, there is no satisfactory 
explanation for the pattern shown in Fig. 3. What is clear, however, is that one 
must be prepared for many surprises when entering the neutron-rich territory. 
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2.2. Isospin Dependence of Pairing 

The uniqueness of drip-line nuclei for studies of effective interactions is due 
to the very special role played by the pairing force. Correlations due to pairing, 
core polarization, and clustering are crucial in weakly bound nuclei. In a drip-line 
system, the pairing interaction and the presence of skin excitations (soft modfislxoiild 
invalidate the picture of a nucleon moving in a single-particle orbit H2P > EJ) > M> > i3 >^3. 

Surprisingly, rather little is known about the basic properties of the pairing 
force. In most nuclear structure calculations, the pairing Hamiltonian has been ap- 
proximated by the state-independent seniority pairing force, or schematic multipole 
pairing interaction©. Such oversimplified forces, usually treated by means of the 
BCS approximation, perform remarkably well when applied to nuclei in the neigh- 
borhood of the stability valley, but they are inappropriate (and formally wrong) 
when extrapolating far from stability The self-consistent mean-field models have 
meanwhile reached such a high level of precision that one needed to improve on 
the pairing part of the model. Presently, the mosLupztQzdate models employ local 
pairing forces parametrized as contact interactions EJ)'0)>LL2P. More flexible forms at- 
tach a density-dependence to the pairing strength©'©'©. There exist even more 
elaborate forms which include also gradient termslHirEj). Although all these various 
density dependencies of pairing are widely used, very little is yet known about their 
relations to observable quantities. 

Up to now, the microscopic theory of the pairing interaction has only seldom 
been applied in realistic calculations for finite nuclei. A "first-principle" derivation 
of pairing interaction from the bare NN force still encounters many problems such 
as, e.g., treatment of core polarization©' 0). Hence, phenomenological density- 
dependent pairing interactions are usually introduced. It is not obvious, how the 
density dependence should be parametrized although nuclear matter calculations 
and experimental data on isotope shifts strongly suggest that pairing is a surface 
phenomenon, and that pairing interaction should be maximal in the surface region. 
This is why neutron-rich nuclei play such an important role in this discussion. Indeed, 
because of strong surface effects, the properties of these nuclei are sensitive to the 
density dependence of pairing. 

Recent workt^iP, based on the spherical Skyrme-HFB model, contains the the- 
oretical analysis of particle and pairing densities in neutron-rich nuclei and their 
dependence on the choice of pairing interaction. In the particle-particle (p-p) chan- 
nel, the density-dependent delta interaction (DDDI)®'©'0 has been employed: 

< P ir (r,r / ) = / P a ir (r)5(r-r'), (2-1) 

where the pairing-strength factor is 

/ P ai r (r) = Vb{l-[pis(r)/p c n (2-2) 

and Vq, p c , and a are constants. The presence of the density dependence in the 
pairing, channel has consequences for the spatial properties of pairing densities and 
fields^'ED'HoD. In Eq. @ pis(r) stands for the isoscalar single-particle density 
/Ois( r ) = /On(r)+p p (r). If p c is chosen such that it is close to the saturation density, 
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Pc~Pis( r = 0), both the resulting pair density and the pairing potential h(r) are 
small in the nuclear interior, and the pairing field becomes surface-peaked. By 
varying the magnitude of the density-dependent term, the transition from volume 
pairing to surface pairing can be probed. 

Apart from rendering the pairing weak in the interior, the specific functional 
dependence on pis used in Eq. (|2-2| ) is not motivated by any compelling theoretical 
arguments or calculations. In particular, values of power a w ere chosen ad hoc to 
be either equal to 1 (based on simplicity), see, e.g., Refs.S'EiP, or^qual to the 
power 7 of the Skyrme- force density dependence in the p-h channelt3)'E9). As a 
typical example, the particle and pairing local HFB+SLy4 neutron densities p n (r) 
and p n { r ) calculated for several values of a are displayed in Fig. 4 for 150 Sn. 
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Fig. 4. Self-consistent spherical HFB+SLy4 local densities p(r) (top) and p(r) (bottom) for neu- 
trons in 150 Sn for several values of a. The insets show the same data in logarithmic scale. 

With decreasing a, the p-p density p n (r) develops a long tail extending towards 
large distances. This is a direct consequence of the attractiveness of DDDI at low 
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densities when a is small. While in the nuclear interior, the p-h density p n (r) depends 
extremely weakly on the actual form of pairing interaction. Due to the self-consistent 
feedback between particle and pairing densities, the asymptotic values of p n { r ) are 
significantly increased when a gets small (see inset). Moreover, one observes a clear 
development of a halo structure, i.e., a smooth exponential decrease, that for a=l 
starts at r~6 fm, is interrupted at r~9 fm for small a, and replaced by a significantly 
slower decrease of the density. The general conclusion drawn from Fig. 4 is that 
experimental studies of neutron distributions in nuclei are extremely important for 
determining the density dependence of pairing interaction in nuclei. 

While the analysis presented in Ref. ® strongly suggested that the strong low- 
density dependence of pairing force, simulated by taking very small values of a in 
DDDI, is unphysical, it is only very recently that a realistic fit of DDDI to the 
odd-even staggering in nuclear masses has been carried outE^P. 
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Fig. 5. Comparison between the experimental staggering parameters (upper left panel, based on 
masses from Ref.B>) and the average neutron pairing gaps calculated within the spherical 
HFB method for the Skyrme SLy4 force and three different versions of the zero-range pairing 
interaction. 

Results of the spherical coordinate-space HFB calculations for semi-magic even- 
even nuclei for the average pairing gaps are shown in Figs. 5 and || for neutrons and 
protons, respectively. In the upper left panels we show the values of experimental 
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three-point staggering parameters Z\® centered at odd particle numbers©'© and 
averaged over the two particle numbers adjacent to the even value. The experimental 
data from the interim 2001 atomic mass evaluation c3 were used. 
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Fig. 6. Same as in Fig. 5 except for the average proton pairing gaps. 

The lower left and right panels in Figs. 5-|6| show the results obtained for the 
surface (a=l) and volume pairing interactions, respectively. When compared with 
the experimental numbers, one sees that both types of pairing interaction are in clear 
disagreement with experiment. The surface interaction gives the pairing gaps that 
increase very rapidly in light nuclei, while the volume force gives the values that 
are almost independent of A. The surface pairing in light nuclei is so strong that 
pairing correlations do not vanish in doubly magic nuclei such as 16 O or 40 Ca. The 
experimental data show the trend that is intermediate between surface and volume; 
hence, one may consider the intermediate-character pairing force that is half way in 
between, i.e., it is defined as: 



i 



P( r ) 
2po 



<5(r-r'). 



(2-3) 



The upper right panels in Figs. 5 and ^ show the results obtained with the mixed 
pairing force. It can be seen that one obtains significantly improved agreement with 
the data, although a more precise determination of the balance between the surface 
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Fig. 7. Similar as in Fig. 5, but for chains of isotones ajjound N=S2. The experimental neutron 
pairing gaps were estimated km the three-pointa'u- 1 staggering parameters centered at 
JV=81 and iV=83 (from Ref.B>). 

and volume contributions still seems to be possible. One should note that similar 
intermediate-character pairing forces have recently been studied in Ref.ES' where it 
was concluded that pairing in heavy nuclei is of a mixed nature. 

Figure 7 illustrates the role of using different types of the pairing interaction 
to predict the neutron pairing gaps in very neutron-rich nuclei. The experimental 
data that exist for Z>50 do not indicate any definite change in the neutron pairing 
intensity with varying proton numbers. However, the surface pairing interactions 
(bottom panels) give a slow dependence for Z>50 that is dramatically accelerated 
after crossing the shell gap at Z=50. On the other hand, the volume and mixed 
pairing forces predict a slow dependence all the way through to very near the neutron 
drip line. It is clear that measurements of only several nuclear masses for Z<50 will 
allow us to strongly discriminate between the pairing interactions that have different 
space and density dependencies. 

§3. Physics of Very Heavy Nuclei 

The stability of the heaviest and superheavy elements has been a longstanding 
fundamental question in nuclear science. Theoretically, the mere existence of the 
heaviest elements with Z>102 is entirely due to quantal shell effects. Indeed, for 
these nuclei the shape of the classical nuclear droplet, governed by surface tension 
and Coulomb repulsion, is unstable to surface distortions driving these nuclei to 
spontaneous fission. That is, if the heaviest nuclei were governed by the classical 
liquid drop model, they would fission immediately from their ground states due to 
the large electric charge. However, in the mid-sixties, with the invention of the shell- 
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correction method, it was realized that long-lived superheavy elements (S HE ) with 
very large atomic numbers could exist due to the strong shell stabilization ElP'Ej). 

The last three years (1999-2001) have brought a number of experimental sur- 
prises which have truly rejuvenated the field (see Ref. for a recent review). Most 
significant are reports fromDuhna^ announcing the discovery of elements 114 and 
116 in hot fusion reactionsl2ZP't^P'E3). Further experiments with stable beams are 
planned; they will be extremely helpful for the theoretical modeling of the SHE for- 
mation. It is anticipated, however, that future experimental progress in the synthesis 
of new elements will be possible - thanks to radioactive nuclear beams, especially 
the doubly magic neutron-rich beam of 132 Sn@. 

In spite of an impressive agreement with experimental data for the heaviest 
elements, theoretical uncertainties are large when extrapolating to unknown regions 
of the nuclear chart. In particular, there is no consensus among theorists with regard 
to the center of the shell stability in the superheavy region. Since in these nuclei the 
single-particle level density is relatively large, small shifts in the position of single- 
particle levels (e.g., due to the Coulomb or spin-orbit interaction) can be crucial 
for determining the shell stability of a nucleus. But does it actually make sense to 
talk about "magic superheavy nuclei"? Recent theoretical workc2P sheds new light 
on this question. According to calculations, the patterns of single-particle levels are 
significantly modified in the superheavy elements. Firstly, the overall level density 
grows with mass number A, as oc A 1 ^ 3 . Secondly, no pronounced and uniquely 
preferred energy gaps appear in the spectrum. This shows that shell closures which 
are to be associated with large gaps in the spectrum are not robust in superheavy 
nuclei. Indeed, the theory predicts that beyond Z = 82 and N = 126 the usual 
localization of shell effects at magic numbers is basically gone. Instead the theory 
predicts fairly wide areas of large shell stabilization without magic gaps (see Fig. 8). 
This is good news for experimentalists: there is a good chance to reach shell-stabilized 
superheavy nuclei using a range of beam-target combinations. 

The Coulomb and nuclear interactions act in opposite ways on the total nucleonic 
density in the nuclear interior. As a consequence of their saturation properties, 
nuclear forces favor values of the internal density close to the saturation density of 
nuclear matter. On the contrary, since the Coulomb interaction tends to increase 
the average distance between protons, the Coulomb energy is significantly lowered by 
either the creation of a central depression or by deformation, or both. Based on this 
general argument, one expects the formation of voids in heavy nuclei in which the 
Coulomb energy is very largeeiP'E3). Recently, the subject of exotic (bubble, toroidal, 
band-like) configurations in nuclei with very large atomic numbers has been revisited 
by self-consistent calculations. The important question which is asked in this context 
is: What are the properties of the heaviest nuclei that can be bound (at least, in 
theory), in spite of the tremendous Coulomb force? 

Calculations do predict the existence of bubble nuclei©-®-© which are stabi- 
lized by shell effects . Figure Jl shows the results of the coordinate-space SLy6+HF 
calculations with delta pairingciP. The potential energy surfaces (PES) correspond- 
ing to various density distributions are plotted as functions of the mass quadrupole 
deformation /?2- This figure nicely illustrates the competition between various forms 
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Fig. 8. Total shell energy (sum of proton and neutron shell corrections) calculated for spherical 
even-even nuclei with 40<Z<60 (top) and for superheavy nuclei around the expected island 
of stability around._Z=120, 7V = 180. Dots mark nuclei predicted to be stable with respect to 
/3-decay (from Ref.EP). 



of nuclear matter constituting a finite heavy nucleus possessing a huge electric charge. 
It is difficult to say at present whether these exotic topologies can occur as metastable 
states or whether they can form isolated islands of nuclear stability stabilized by shell 
effects. The interplay between exotic configurations in superheavy nuclei is likely to 
impact the exact position of the borders of the superheavy and hyperheavy region. 

§4. Conclusions 

The main objective of this brief review was to discuss various facets of RNB 
physics. The list of topics covered is, of course, very limited due to time and space 
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Fig. 9. Variety of shapes predicted for the superheavy nucleus 780 2 54526 in the SLy6+HF model 
with delta pairing. Axial symmetry was assumed. The local minimum at /?2=0 can be associated 
with a bubble configuration (solid line). At small deformations, the PES corresponding to the 
"normal" density distribution (dashed line) is predicted to lie very high in energy. At small 
oblate deformations, the unstable "band" -like structure becomes favored energetically. The 
single-particle levels corresponding to the spherical bubble minimum are shown in the inset. 
The contour plots of the corresponding total densities are given in the boxes. The contours are 
axially symmetric with reeeect to r=0 and the contour lines are at 0.01, 0.03, 0.06, 0.09, 0.12, 
and 0.15 fm" 3 (from Ref.EP). 

constraints. An experimental excursion into uncharted territories of the chart of 
the nuclides, exploring new combinations of Z and N, will offer many excellent 
opportunities for nuclear structure research. What is most exciting, however, is that 
there are many unique features of exotic nuclei that give prospects for entirely new 
phenomena likely to be different from anything we have observed to date. We are 
only at the beginning of this most exciting journey. 
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